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ABSTRACT 

The propagation of dissipationless, hydromagnetic, purely toroidal Alfven waves in a realistic back- 
ground three-fluid solar wind with axial symmetry and differential proton-alpha flow is investigated. 
The short wavelength WKB approximation is not invoked. Instead, the equations that govern the 
wave transport are derived from standard multi-fluid equations in the five-moment approximation. 
The Alfvenic point, where the combined poloidal Alfven Mach number Mt = 1, is found to be a 
singular point for the wave equation, which is then numerically solved for three representative angular 
frequencies u> = 10~^, 10"'' and 10~^ rad s~^ with a fixed wave amplitude of 10 km s~^ imposed at the 
coronal base (1 Rq). The wave energy and energy flux densities as well as wave-induced ion accelera- 
tion are computed and compared with those derived in the WKB limit. Between 1 Rq and 1 AU, the 
numerical solutions show substantial deviation from the WKB expectations. Even for the relatively 
high frequency lo = 10^^ rad s~^, a WKB-like behavior can be seen only in regions r > 10 Rq. In the 
low- frequency case w = 10~^ rad s~^, the computed profiles of wave-related parameters show a spatial 
dependence distinct from the WKB one, the deviation being particularly pronounced in interplanetary 
space. In the inner corona r < 4 Rq, , the computed ion velocity fluctuations are considerably smaller 
than the WKB expectations in all cases, as is the computed wave-induced acceleration exerted on 
protons or alpha particles. As for the wave energy and energy flux densities, they can be enhanced 
or depleted compared with the WKB results, depending on oj. With the chosen base wave amplitude, 
the wave acceleration has negligible effect on the ion force balance in the corona. Hence processes 
other than the non-WKB wave acceleration are needed to accelerate the ions out of the gravitational 
potential well of the Sun. However, at large distances beyond the Alfvenic point, the low-frequency 
waves can play an important role in the ion dynamics, with the net effect being to equalize the speeds 
of the two ion species considered. 

Subject headings: waves — Sun: magnetic fields-solar wind-Stars: winds, outfiows 



1. INTRODUCTION 



Eve r since their identification by iBelcher fc David 
(flQTlh . Alfven waves have been extensively studied 
using in situ measurements, such as by Helios and 
Ulysses, cov ering the hel i ocent r ic distance from 0.29 
to 4.3 AU (iTu fc MarschI [l995l : iGoldstein et"an fl995l : 
iBavassano et al.ll2000a rbl). On the other hand, the non- 
thermal broadening of a number of Ultraviolet lines, such 
as those measured with the SUMER (Solar Ultraviolet 
Measurements of Emitted Radiation) and UVCS (Ultra- 
violet Coronagraph Spectrometer) instruments on SOHO 
(the Solar and Heliospheric Observatory), is usually at- 
tributed to the transverse velocity fluctuations, thereby 
enabling one to infer the amplitudes of these fluctuations 
in the inner corona below ~ 5 Rq (Bancriec et al. 1998; 
lEsser et al.lll999( ). Moreover, the Faraday rotation mea- 
surements, which yield information regarding the line- 
of-sight magnetic field fluctuations, have been shown to 
support indirectly the presence of Alfven waves inside 
10 Rq (HoUweg et al. 1982). The solar wind in interme- 
diate regions, for the time being, can be explored only 
by radio scintillation measurements which allow one to 
derive the velocit y fluc tuations (e.g.. [Armstrong fc Wool 
119811: IScott et al.lll983D . It is noteworthy that although 
the hourly-scale fluctuations seem to be more frequently 
studied, the fluctuation spectrum measured by Helios 
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nevertheless spans a broad frequency range from 10 ^ 
to 10-2 (lTu_&JvIarsch 1995). 

Most of the theoretical investigations into the interac- 
tion between Alfven waves and the solar wind have been 
performed in the short wavelength WKB limit which 
makes the problem more tractable mathematical ly. For 
instan ce, by employing the WKB approximation. IParkeil 
(|1965D derived an expression for the pondcromotive force 
through which the Alfven waves may provide further ac- 
celeration to the solar wind. The wave acceleration was 
later incorporated in detailed numerical models by, e.g., 
Alazraki fc Couturier (1971). It was soon realized that 
Alfven waves may also heat the solar wind via dissipa- 
tive processes such as the cyclotron resonance interac- 
tion between ions and high frequency, parallel propa- 
gating waves genera ted by a turbulent c a scade (cf. the 
extensive review bv iHollw eg fc Isenberd (|2002t )). Such 
a parallel cascade scenario has been successful in ex- 
plaining a number of observations, to name but one, the 
significant thermal ani sotropy of ions as established by 
U yCS measurements (ITet al.lll999| ). As pointed out 
by IHollweg fc Isenberg ( 2002f ). the applicability of the 
WKB approximation in which the processes are formu- 
lated is questionable in the near-Sun region in view of the 
large Alfven speeds. Furthermore, a turbulent cascade 
requires a non-vanishing magnetic Reynolds stress ten- 
sor, which, however, is zero in the WKB limit since the 
particle and field components of the tensor cancel each 
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other exactly. A non-WKB analysis is therefore required 
to account for the wave reflection and the consequent 
driving of any turbulence cascade. 

As a matter of fact, non-WKB analysis of Alfven waves 
in the solar wind has been carri ed out for decades (e.g., 
iHeinemann fc Olbertlll980l : lLoullT99aV This, however, is 
almost exclusively done in the framework of ideal MHD, 
which allows waves propagating in opposite directions to 
be explicitly separated when the Elsasser variables are 
used. The adoption of the Elsasser variables has also en- 
abled a new turbulence phenomenology concerning the 
nonlinear coupling b etween counter-propagating waves 
(iDmitruk et al.l[2n0ll :[C ranmer fc van Ballegooij enl [20?)5l : 
iVerdini et al.ll2005D . This coupling term, if interpreted 
as the energy cascaded towards fluctuations with increas- 
ingly large perpendicular wavenumbers, is also more con- 
sistent with theoretical expectations. 

Despite substantial advances achieved in ideal MHD, 
which is appropriate for the description of the gross 
properties of the solar wind, the non-WKB analysis 
of Alfven waves has rarely been done using multi-fluid 
transport equations. In contrast, multi-fluid, Alfven 
wave driven solar wind modeling formulated in the 
WKB limit has reached c onsiderable sophistication (see 
iHollweg fc Isenberell2002l ). Such a multi- fluid approach 
is particularly necessary for the solar wind since the al- 
pha particles must be included given their non-negligible 
abundance and the fact that the proton-alpha differential 
speed can be a s ubstantial fraction o f the proton speed in 
the fast stream (jMarsch et al.l[l982l ). There is therefore 
an obvious need to extend available non-WKB analyses 
of Alfven waves from the ideal MHD to the multi-fluid 
case. 

The intent of this paper is to present an analysis of 
Alfven waves in a 3-fluid solar wind assuming axial sym- 
metry, without assuming that the wavelength is small 
compared with the spatial scales at which the back- 
ground parameters vary. The perturbed velocity and 
magnetic field are assumed to be in the azimuthal direc- 
tion, i.e, only purely toroidal waves are investigated. To 
further simplify the treatment, the wave dissipation is ne- 
glected. This simplification is necessary here since if one 
wants to gain some quantitative insights into the wave 
dissipation, and to maintain a reasonable self-consistency 
at the same time, one should formulate the dissipation in 
terms of the amplitudes of waves propagating outwards 
and inwards. To this end, a full multi-fluid Elsasser anal- 
ysis is required but is unfortunately unavailable at the 
present time. However, if the wave dissipation is ne- 
glected and therefore the nonlinear interaction between 
waves and the multicomponcnt wind is entirely through 
the agent of ponderomotive forces, then the problem can 
be formulated without distinguishing explicitly b etween 
the directions of wave propagation fcf. lLouiri993[ ). 

The paper is organized as follows. In section[51 we show 
how to reduce the general multi-fluid transport equations 
to the desired form. The resultant equations govern- 
ing the Fourier amplitudes of toroidal Alfven waves at a 
given frequency are then solved analytically in two limit- 
ing cases, namely the WKB and zero-frequency limits, in 
section [3l Apart from these analytically tractable cases, 
the equations have to be solved numerically. In sectionlH 
we re-formulate the model equations for numerical con- 
venience, describe the background flow parameters and 



detail the solution procedure as well. The numerical so- 
lutions for three different frequencies are presented in 
section [SJ Finally, section [S] summarizes the results, end- 
ing with some concluding remarks. 

2. MATHEMATICAL FORMULATION 

Presented in this section is the mathematical develop- 
ment of the equations that govern the toroidal fluctu- 
ations in a solar wind which consists of electrons (e), 
protons {p) and alpha particles (a). Each species s 
(s = e,p,a) is characterized by its mass rris, electric 
charge e^, number density n^, mass density ps = nsUis, 
velocity Vj, and partial pressure Ps- If measured in units 
of the electron charge e, Cg may be expressed by — ZgC 
with Zg = — 1 by definition. 

To simplify the mathematical treatment, a number of 
assumptions have been made and are collected as follows: 

1. It is assumed that the solar wind can be described 
by the standard transport equations in the five- 
moment approximation. 

2. Quasi neutrality and quasi-zero current are as- 
sumed, i.e.. Tig = J2k^knk and Vg = 
J2k Zknk^k/ne where k ^ p,a. 

3. Symmetry about the magnetic axis is assumed, i.e., 
d/d4> = in a heliocentric spherical coordinate sys- 
tem (r, 0, (j)). 

4. The time-independent solar wind interacts with the 
waves only through the wave-induced ponderomo- 
tive forces. 

5. The wave frequency considered is in the hydro- 
magnetic regime, i.e., well below the ion gyro- 
frequencies. 

6. The perturbed velocity and magnetic field are as- 
sumed to be in the (p direction only. 

7. The effects of the solar rotation on the background 
solar wind are neglected such that there is no need 
to consider the coupling of Alfven waves to the 
compressional modes in the presence of a spiral 
magnetic field. 

8. The effects of the Coulomb friction on the waves are 
neglected, and so is the wave-induced modification 
of the Coulomb friction between background ion 
flows. 

2.1. Multi-fluid Equations 

The equations appropriate for a multi-component solar 
wind plasma in the standard five-moment approximation 
are as f ollows (for the derivation see appendix A.l in 
iLi fc Lil (11)06)) 

V • (nfeVfe) 



dt 
dt 



Vfc • Vvfc 
1 



= 0, 
Vpfc 



(1) 
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ZkG TljZj 



V,, • V- 



(vj - Vfc) X B = 0, 



7 — 1 7^1 
SE,, 



+ V • q, - ^ - = 0, 

— - V X (Ve X B) = 0, 



(2) 

(3) 

(4) 



where the subscript s refers to ah species (s — e,p,a), 
while k stands for ion species only {k = p,a). The grav- 
itational constant is denoted by G, is the mass of 
the Sun, B the magnetic field and c the speed of light. 
The momentum and energy exchange rates due to the 
Coulomb collisions of species s with the remaining ones 
are denoted by S'M.s/St and SEg/St, respectively. More- 
over, Qs is the heat flux carried by species s, and Qs 
stands for the heating rate applied to species s from 
non-thermal processes. In equation the subscript 
j stands for ion species other than fc, namely, j — p for 
k = a and vice versa. As can be seen, in addition to 
the term (V x B) x B, the Lorentz force possesses a new 
term in the form of the cross product of the ion velocity 
difference and magnetic field. Physically, this new term 
represents the mutual gyration of one ion species about 
the other, the axis of gyration being in the direction of 
the instantaneous magnetic field. 

Equations ([T|) to ([3]) form a complete set if supple- 
mented with the description of species heat fluxes q^ 
and heating rates Qs- As such, they can be invoked 
to depict self-consistently the interaction between Alfven 
waves and the solar wind species by explicitly introducing 
these waves via boundary conditions. On the one hand 
solving this set of equations presents a computationally 
formidable task; on the other hand, one can extract the 
necessary information concerning the dynamical feed- 
back of the waves to the plasma by going beyond the 
WKB limit. In the non-WKB approach to be adopted 
here, one assumes that the governing equations can still 
be separated into those governing the background time- 
independent flow a nd those g overning the transport of 
waves. As noted bv lLoul (|1993l ) (also see the discussion), 
this separation does not necessarily require the waves be 
linear as long as sufficiently small wave amplitudes are 
imposed at the Sun. 

Further simplification also results from the choice of a 
flux tube coordinate system, in which the base vectors 
are {e;, bat, e^}, where 



Bp/Bf 



ejv = 60 X ei, 



with the subscript P denoting the poloidal component. 
Moreover, the independent variable I is the arclength 
along the poloidal magnetic field line. This choice per- 
mits the decomposition of the magnetic field and species 
velocities into background ones and fluctuations. 



B = Biii + he. 



= UsCi + WsCn + UsB^,. (5) 



where s = e,p,a. From the assumption of az- 
imuthal symmetry, and the assumption that Bp is time- 
independent, one can see from the poloidal component 
of equation (|3|) that v^p should be strictly in the di- 
rection of Bp. In other words, We = to a good 
approximation. Now let us consider the component 



of the momentum equation ([2]). Since the wave fre- 
quencies in question as well as other frequencies asso- 
ciated with the spatial dependence are well below the 
ion gyro-frequency fifc = {ZkcBi) / {mkc) (k = p,a), 
from an order-of-magnitude estimate one can see that 
\wj — Wk\ <C \uk\- Combined with the fact that We = 0, 
this leads to that both Wp and Wa should be very small 
and can be safely neglected unless they appear alongside 
the ion gyro-frequency. With this in mind, one can find 
from the N component of equation ^ that 



IT (f/o - Up) . 



(6) 



That is, the ion velocity difference is aligned with the 
instantaneous magnetic field. This alignment condition 
further couples one ion species to the other. Note that 
due to the assumption of quasi-zero current, equation ([6]) 
leads to 



Uk 



(7) 



where k — p,a. 

Given the aforementioned assumptions, the time- 
independent multicomponent solar wind in which the 
toroidal Alfven waves propagate is governed by 



n-kUk 
B, 



= 0, 
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6t 



(Inr)' 
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SE, 
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-B, 



(8) 



(9) 



(10) 



where the prime denotes the derivative with respect to 
the arclength I which becomes the only independent spa- 
tial variable. In addition, a^ k denotes the acceleration 
exerted on ion species k {k = p, a) by the toroidal fluc- 
tuations. 



aw,k = (ul) (Ini?)' - 
{bXk) 



ATmemk 



B, 



(11) 



where R — r sin 6 is the distance from a point along the 
poloidal magnetic field line to the magnetic axis. Besides, 
the angular brackets stand for the time-average over one 
wave period. The variable Xk — ^k{wj — Wk){Zjnj /ue) 
distinguishes the present study fror n those using the ideal 
MHD in which case Xk = (e.g., iHeinemann &: Olbertj 
Il980f ). Apart from this, the wave-induced accelera- 
tion a^^k includes the inertial centrifugal acceleration 
(the first term on the right hand side), and the usual 
(V X B) X B term (the second term). 

2.2. Transport of Toroidal Alfven Waves 

The transport of toroidal Alfven waves is governed by 
the azimuthal component of the momentum equation ([2]) 
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for ion species k [k — p, a) together with that of the 
magnetic induction law To be more specific, these 
equations read 



1 fUe 



duk 
dt 



Uk 



duk 
dl 



R \Bi 
Uk{\nR)' 



0, 



Zk 



-B, 



db 



5(lni?)' 



Xk. 



(12) 



(13) 



The wave propagation is characterized by several key 
parameters, namely, the wave energy and energy flux 
densities as well as the wave-induced acceleration. These 
parameters can be found by considering the energy con- 
servation for the Alfven waves, 



2 Stt^ ^dl\Bi 



E 

L k 



pkUkul 



An 



(Ueb - BiUe 



Ue 



Y^PkUkulilnR)' , 

fc 

y^PkUkXk. 



b—b + b^ilnRY 
ol 



(14) 



On the left hand side (LHS), the first term is the time 
derivative of the perturbation energy density, while the 
second is the divergence of the perturbation flux density. 
The physical meaning of the right hand side (RHS) can 
be revealed as follows. By using relation ([7]) one finds 



PkUkXk - ^ PkUk 



bXk 
Bi 



Bi 



Since J2k PkXk = 0, one can identify the time-average 
of the RHS of equation (HH) as the negative of the work 
done by the wave- induced forces on the solar wind, i.e., 
J2k PkUkaw,k (cf. Eq. lfTTjl V In other words, taking the 
time-average of equation (|14p yields 



B4l^ = -T.P'^Uka^ 

' A: 



w.k 1 



(15) 



where is the time-average of the perturbation flux 
density and will be termed wave flux density for simplic- 
ity. (Actually, the wave properties to be discussed always 
refer to time-averaged values.) The total energy is there- 
fore conserved for the system comprised of a multi-fluid 
solar wind and toroidal Alfven waves: The gain in the 
solar wind kinetic energies is at the expense of the wave 
energy. 

The appearance of Xk makes equation (|13p inconve- 
nient to work with. Instead one may consider the az- 
imuthal component of the total momentum to eliminate 
Xk, the resultant equation being 



k 

B. 



d 



AttR dl 



duk Uk d ■ 



[Rb) = 0. 



(16) 



Now one may proceed by introducing the Fourier ampli- 
tudes at a given angular frequency w, 

[5(/, t),uS: t)] = W). us{l)] exp(-ic^i), (17) 

in which s — e,p, a. As a result, the Faraday's law ([1 
and equation (jl6[) now take the form 



Uj)' - Biu' 



iuj + Ue In 



RBi 

Ue 



b- Bi{\nR)'{Le 



(18) 



^T.PkUk'^-^Yb'^Y.PM'e 

\ k ' Ik 

^^('i?^^^V-^(lni?)' 



R 



B, 



Y^PkUkihiR)''. 



. ^ ( Uk-Ue l , , 

ibj^Pk I — ^ -I- u. 



(19) 



in which we have used relation ^ to express the ion 
velocity fluctuation Uk in terms of the electron one Ue- 

Once the background flow parameters and proper 
boundary conditions are given, equations p8)) and p9|) 
can be solved for the Fourier amplitudes of magnetic fluc- 
tuation b and the electron velocity fluctuation Ue for a 
given angular frequency uj. The ion velocity fluctuations 
Uk can then be found in virtue of relation (Q). The three 
wave-related parameters can be evaluated by forming a 
time-average over one wave period. Speciflcally, the wave 
energy density i^u, and wave energy flux density F^^, are 
given by 



E 

k 

E 



Pk{ul) . {b') 



(20a) 



PkUk {u\) 



— {Ue (b^) - Bi {bUe)) 



(20b) 



while the wave-induced acceleration aw,k has already 
been given by equation (fTTj) . It can be seen that 
consists of the ion kinetic as well as the magnetic ener- 
gies, while F^ is comprised of the kinetic energy flux con- 
vected by ion fluids and the Poynting flux. To evaluate 
the time-average (fg) for two wave-related fluctuations 
f.fi.gl f.fi.9] exp(— ia;^), one may use the expression (cf. 
iMacGregor fc CharbonneaullTOM ) 

(/g) = iRe(/r) = J(/r + f5), (21) 

where Re(z) denotes the real part of a complex variable 
z, and the superscript asterisk denotes the complex con- 
jugate. 

Although in general one has to rely on a numerical 
integrator to solve equations (fT8|) and (fT9|) . as shown 
below, they are analytically tractable for two limiting 
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3. ANALYTICAL SOLUTIONS IN TWO LIMITS 

In this section, we will examine the analytical solu- 
tions to equations ([18]) and (|19p in the WKB and zero- 
frequency {lo = 0) limits. These analytical treatments 
not only help to validate numerical solutions, but also 
allow one to gain insights into the mathematical proper- 
ties of the governing equations. 

3.1. the WKB Limit 

Extensive studies have been made on the Alfven waves 
in the WKB limit. In particular, the Alfven wave force 
exerted on differentially streaming ioi is was first de- 
rived by iHollwed (|1974l ) and later by iMcKenzie et all 
(|1979f ) for a spherically expanding solar wind. Based 
on the idea of the wave av erage Lagrangian, the work 
bv'Isenb erg fc Hollwed (1982) not only obtained the ex- 
pression for the wave force but further showed that there 
exists an adiabatic invariant, namely the wave action 
flux, in the absence of wave dissipation. The deriva- 
tion is rather general and does not require a specific 
flux tube geometry. Alternatively, the actio n flux con- 
servat ion has been independently obtained by iMcKenzi^ 
(|1994[ ) who used the mo re familiar WKB analysis. Al- 
though iMcKenzii (|1994D assumed that the solar wind 
is again spherically symmetric, his results can be read- 
ily extended to an axisymmetric conflguration. In what 
follows, such an extension is presented. 

The formal development of the WKB analysis starts 
with the introduction of the expansion 



Us = 



hUs,2 + ---)exp(iS'(/)), (22a) 
b2 + ---)exp{iS{l)), (22b) 



in which s — e,p,a. It is assumed that Ug^n and &„ 
{n= 1,2,...) vary at the same spatial scale H as the back- 
ground flow parameters, and 1/7? is small compared with 
the wave number K = S'{1), i.e., /i — 1/{KH) <C 1. In 
addition, it is assumed that \us^n+i/us^n\ ^ \bn+i/bn\ ^ 
/i. Substituting the expansion ((22l) into equations (|18p 
and (|19p . one can flnd by sorting different terms accord- 
ing to orders of ^ that 



LJ ehn + KBiUt 

-.KB, ' 



R\B, 



bn-l 



Ue.n— 1 



5/ 



(23a) 



k 



i— 



R 



PkUk 
R 

{Rbn-lY : 



Uk-Ue 



^e,n — 1 



bn 

Uk - u. 



4tt 



B, 



■bn- 



: uj - KUs {s 



(23b) 

e,p, a) has been 



where the definition lJs - 
used. 

At the lowest order n = 1, the RHS of equation (|23p 
is zero. For Me,i and bi not to be identically zero, the 
determinant of the coefficient matrix has to be zero. As 
a result, one finds the local dispersion relation 



(24) 



where Uph = uj/K is the local phase speed of Alfven 
waves. Furthermore, one can find the local eigen- 
relations at order n = 1, 



Us,l = - {Uph - Us) 



h 

Br 



(25) 



with s — e,p, a. 

At an arbitrary order n > 2, it follows from equa- 
tion ((23a|) that 



{Uph -Ue)-f^ + Un, 



in which 



iR 



1 fUe 



R\B, 



arises from quantities at order n — 1. At order n + 1, 
one can eliminate bn+i and Ue^n+i on the LHS of equa- 
tion by using the dispersion relation ([M]) , the resul- 
tant equations taking the form 



Tn + Sn — 0, 



(26) 



where 



Pk^^kRBi f ^ 



PkUk 
R 



Tn — ^ Pkl^kRBl 



{RUn)' 
Uphbn 

RBi 



(27a) 



PkUk ( j^ Uph - Uk , 



R 



Bi 



(27b) 



The barred and unbarred parts are again distinguished. 

On expanding T„, by repeatedly using the dispersion 
relation (|24p one finds that 



Tn = 



Uphbn 



Uphbn\ p{Uph-Um) 



B, 



B, 



(28) 



in which the fo llowing definitions have been used (cf. 
lMcKenzid[l99l 



P = '^Pk, 

k 

pUm = ^ PkUk, 



pU^^Y.P'^^' 



(29a) 
(29b) 
(29c) 



with the last one defined for future use. If putting n = 1, 
from 5*1 — one finds Ti = 0, i.e.. 



^, p{Uph-Un.)U^, 

Oi — ^5 — = const. 

Bf 



(30) 



Note that equation (|30p can be interpreted in terms of the 
conservation of wave action flux, which was first found 
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by llsenberg fc HoUwegj ()1982f ). On the other hand, ex- 
panding equation (fTS]) to first order in ^ and using the 
resultant Xk to evaluate {bXk) in equation pT|) . with the 
aid of the eigen-relation (I25p . one can find a compact ex- 
pression for the wave-induced acceleration on ion species 
k {k = p, a), 



r/2 



2Bf 



(31) 



where the time-average (6^) = |6ip/2. It is noteworthy 
that the shape of the line of force, which enters into the 
discussion through R — rsin^, does not show up via 
Sn-i unless one examines the evolution of higher order 
fluctuations Us^n and 6„ {n> 2). 

Some remarks on the wave acceleration a^^^fc given 
by Equation (j3ip are necessary. One may find that 
{U^ - Ul){l + l&il V(2B2)) = const if the flow speeds are 
entirely determined by the wave forces. Usually \bi\/Bi 
increases with increasing distance. The net effect of au],k 
is therefore to limit the proton-alpha differential speed 
Uap- In this sense flm.fc was liken ed to an additi onal fric- 
tional force (e.g.. iHollwedllQTif ). As noted bv iHoUwe^ 
(|1974[ ). the differential nature of ayj^k derives eventu- 
ally from a wave-induced transverse drift velocity 5ufc 
{k = p,a), which is different for different species be- 
cause they see different wave electric fields when they 
flow differentially. The contribution to the wave force 
from (5ufe is twofold. First, it induces a centrifugal force 
that is always non-negative. Second, it contributes to 
the species electric current flow 5Jj_,fc which in turn ex- 
erts on ion species k a 53±^k x (5B force, where (5B is 
the wave magnetic field. This force may become nega- 
tive. As a matter of fact, the two terms may even add 
up to a negative value: The wave force may decelerate 
rather than accelerate the ion species. Conside r a sim- 
plified situation similar to that considered by HoUwegl 
(jl974( ). One assumes that the spatial variation of both 
Ua and Up can be neglected, and Ua > Up > Ua where 
Ua — Bi/y/Airp is the bulk Alfven speed. One may fur- 
ther assume that the alpha particles are test particles, 
the background magnetic field is purely radial, and that 
(X r"*^ with e being a positive constant. In this case, 
from Equation ([3T|l one finds that 



\bi\ 



ABfr 



{ [(4 - e) C/2 + eUl] - (4 - e) {U^p + Upf] 



It follows that a^u.ct < when Uap > Uc where the critical 
value Uc can be roughly approximated by 



Ua 



Ua 



2(4 



Ua, 



which may be smaller than Ua by a factor of Ua/Up- 
Noting that now the wave phase speed is Uph = Up + Ua, 
one finds that the waves may decelerate alpha particles 
even though Ua < Uph- It turns out that this much 
simplified picture roughly represents the behavior of aw,a 
in the region between 20 and 30 i?© where e can be taken 
to be 2.5 for the adopted background flow parameters 
(see Figure[SjD which indicates that ^ in the WKB case 
becomes negative beyond 22.7 Rq where Uap exceeds 

f/c). 



3.2. the Zero-frequency Limit 
The other extreme will be a> = 0. In this case, equa- 



tions p2| and (fT6|) can be integrated to yield 
b 



Ue 



R 



E 



Bi 

PkU k 

pUr, 



AnR, 



-Uk 



Bib 
AnpUr, 



A, 



(32a) 



(32b) 



where Aq and A^ are two integration constants. When 
deriving equation (|32bp . we have used the fact that 
pUm/Bi is a constant. Equation ([5^ demonstrates a 
clear connection to the problem of angular momentum 
transport in a multicompo nent solar wind, as has been 
discussed bv lLi fc Lil (|2006l) . 

With the aid of equation ([7]), equation (|32p can be 
solved to yield 



R 



AnpUn, Al/R^ - Ai 



where 



B, 



iirpUf 



(33) 



(34) 



is the square of the combined Alfven Mach number. For 
typical solar winds, there exists one point where Mt = 1 
between 1 Rq and 1 AU. Call this point the Alfvenic 
point, and let it be denoted by subscript A. It then 
follows that for b not to be singular at this point, one 
must require A^ = AqR'^^. As a result, one may obtain 



b =AnR 

Us ^AfiR 
where s = e,p, a. 



AirpU^ {Ra/RY - 1 



Bi Ml - 1 ' 

, ATTpU,nUs{RA/Rf -I 



Ml ~ 1 



(35a) 
(35b) 



4. NUMERICAL MODEL AND METHOD OF SOLUTION 

Apart from the WKB and zero- frequency limits, equa- 
tions (jlSp and (fT9|) can be integrated only numerically. 
Nevertheless, the treatment in the zero-frequency limit 
in the previous section reveals that the Alfvenic point 
where = 1 is a singular point for the system of equa- 
tions. One may expect that this is also the case for an 
arbitrary finite w. In their present form, however, equa- 
tions (fT5|) and (fTO|) do not show explicitly the existence 
of such a singular point. Hence they are not convenient 
to work with numerically and need to be further devel- 
oped in what follows. Moreover, we will also describe 
how to specify the background flow parameters, and the 
solution procedure as well. 

4.1. Further Development of Governing Equations 

For the convenience of presentation, one may define, in 
addition to Mr, the following dimensionless parameters 

AirpUmUc ^ AirpUeUe 



_ AnpUeUA Y _ A-KpUmUA 
^eA — ,^mA — -g2. J 



z 



Ml 



Bf 

-i-x„ 



(36) 
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in which Ua — Bi/^Anp is the bulk Alfven speed. Now 
equation (fT9|) can be expressed as 



^r, , 4^PC/™., (RZy- 4^pf/,„(lni?)', 

■^0 H U„ = Up, 



R 



B, 



(37) 



Taking into account equation (jlSp . one finally arrives at 

{M^ ^ I) = Fn^ + F,2V, (38a) 

(M| - 1) ry' = ^^21^ + F22V, (38b) 

where two dimensionless variables have been introduced 
for convenience, 

i = b/Bi, r] = Ue/UA- 
Moreover, the coefficients are 

R 

fIX^„ - A„„) + A„,„ I In 

Ue 



(39) 



i^ii 



{RBizy 

RBi 



F, 



12 



U, 



■XeA-2X^A (Ini?)' 



F21 — —j-r- [Z — X„ip + Xe 



F22 — 



Ua 
-Xe 

Ue 

(M| - 



Z In 



me 

RBi\' , {Rzy 



Ue 



R 



+ {Z-X„,e){\nRy 
l)(lnC/A)'. 



(40a) 
(40b) 

(40c) 
(40d) 



It is now clear that the Alfvenic point where = 1 
is a singular point for equation psp . The ensuing task 
is therefore to find solutions that pass smoothly through 
the Alfvenic point. 

4.2. Background Three-fluid Solar Wind 

In principle, one needs to solve equations dH) to (fTO| 
to find a realistic background solar wind by neglecting 
the wave-induced acceleration in equation ^ and us- 
ing a suitable heat input Qs (s = e,p, a). However, 
it is observationally established that in the heliocentric 
range r > 0.3 AU, the proton-alpha differential speed 
Uap = Ua — Up closely tracks the local Alfven speed in 
the fast solar wind with Up > 600 km s~^ (|Marsch et al] 
|1982j). So far this fact still poses a theoretical challenge: 
Adjusting the heating parameters proves difficult to yield 
such a behavior for Uap. In what follows, we shall adopt 
a two-step approach to find the background flow param- 
eters. First, equations ^ to PH)) are solved by using 
a suitable set of heating parameters along a prescribed 
meridional magnetic field line. As a result, we have pp. 
Up, {paUa)i and Uap.i for the whole heliocentric range 
from 1 Rq out to 1 AU. Note that the subscript / de- 
notes the results obtained from this step. Second, the 
alpha parameters are altered as follows. We first impose 
an ad hoc profile for Uap which is identical to Uapj in 
the range r < 0.3 AU but undergoes a smooth transi- 
tion to a profile that is roughly 0.8 Ua everywhere for 




L 






P 

(b) - 




Fig. 1. — Background 3-fiuid solar wind in which the 
toroidal Alfven waves propagate, (a) Adopted poloidal magnetic 
field configuration extending from the coronal base (1 -Rq) to 
1 AU=215 Rq. Note that the magnetic axis points upwards. The 
thick contour delineates the line of force along which the wave 
equation is integrated, (b) Proton and alpha flow speeds Up and 
Ua (solid lines), together with the phase speed Up^ (dashed line) 
expected in the WKB limit as given by equation I l24t . (c) Proton- 
alpha differential speed Uap = Ua — Up (solid line) and the lo- 
cal bulk Alfven speed Ua = -Bj/V^tP (dashed), where Bi is the 
poloidal magnetic field strength and p the overall mass density. 
The asterisks in panels (b) and (c) refer to the Alfvenic point where 
Afy = 1, My being the combined Alfven Mach number defined by 
equation 1134 l l. 



r > 0.3 AU. The alpha speed and mass density are then 
given by Ua = Up + Uap and pa = {paUa)i /Ua, respec- 
tively. 

For the meridional magnetic field, we will ad opt an 
analytical model given by ,Banaszkiewicz et al.l (fl998,) . 
In the present implementation, the model magnetic field 
consists of dipole and current-sheet components only. A 
set of parameters M ^ 3.6222, Q = 0, A' = 1.0534 and 
ai = 2.5 are chosen such that the last open magnetic 
field line is anchored at heliocentric colatitude 6 = 40° 
on the Sun. 

The background magnetic field configuration and flow 
parameters are depicted in Figure [TJ In FiglT^, the field 
line along which we integrate the solar wind equations ([5]) 
to ([To]) and the wave equation ([38)1 is delineated by the 
thick contour. This line of force is rooted at colatitude 
9 = 31.5° on the Sun where the poloidal magnetic field 
strength Bi is 6.6 G. It reaches = 70° at 1 AU where 
Bi is 3.37, compatible with the Ulysses measurements 
(ISmith fc Baloghlll995l ). In FigHb, the sohd lines give 
the distribution with heliocentric distance r of the proton 
[Up) and alpha [Ua) speeds, while the dashed curve is for 
the phase speed [Uph) of Alfven waves in the WKB limit 
as determined by equation ([M)) . Furthermore, Fig[T}: 
shows the distribution with r of the proton-alpha differ- 
ential speed Uap (the solid line) and Alfven speed Ua 
(the dashed curve). The asterisks in FigsUb and [It de- 
note the Alfvenic point where M|, — 1, which is located 
at rA = 11.6 Rq in the chosen background model. 

It can be seen in Fig[T)D that Upt is much larger than 
the ion flow speeds below ~ 20 Rq but is close to Ua 
beyond ~ 0.3 AU. This latter behavior is understandable 
since one can find from equation ([24| that 



Uph Uj-f 



PaPp jj2 
r,2 



1 Xrj 



(41) 
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where Xk = {pkl p){Uap/UA) {k = p,a). The second 
expression is very accurate for Uap < Ua which is the 
case for the assumed background flow (cf . FiglT]:) . Hence 
Uph ~ Ua in the region r > 0.3 AU where Ua > Ua and 
Xp + XpXa/2 ~ 1. Moreover, it can be seen from Fig[Tt 
that Uap drops from nearly zero at Rq to —38.7 km s~^ 
at 1.42 Rq, beyond which Uap increases graduafly. In 
interplanetary space r > 0.3 AU, Uap can be seen to 
follow Ua closely, as is required. 

At the Alfvenic point va, the proton speed Up is 
602 km s~^, the alpha one Ua is 721 km s~^. Eventually 
Up reaches 648 km s^^ at 1 AU where Ua is 676 km s^^. 
As for the density parameters rip and Ha, the model 
yields Up = 3.57 cm~'^ and Ua/up — 4.6% at 1 AU. 
As a result, the proton (alpha) flux UpUp (uaUa) is 2.32 
(O.ll)xlO^ cm"^ s"^ when scaled to 1 AU. Afl these 
values are consistent with typical measurements for the 
low-latit ude fast solar wind str eams, e.g., those made by 
Ulysses (jMcComas et al.ll2000D . 

4.3. Solution Procedure and Boundary Conditions 

Given the background flow parameters, the coefficients 
in equation can be readily evaluated for a given an- 
gular frequency ui. Equation ([38|) is ready to solve once 
proper boundary conditions are supplemented. As is well 
known, one boundary condition has to be imposed at the 
Alfvenic point to ensure the solution is regular. To es- 
tablish this, let us consider the region adjacent to the 
Alfvenic point Ia, such that any function f{l) can be 
Taylor-expanded , 

/(o = + a -^a), (42) 

in which /^"^ — /{Ia)- When substituting this expansion 
into equation p8p . one finds 

F(?e(°)+^^(°)r;(«)=0, (43a) 
F^f^('^+F^%('^^0, (43b) 
F^o)^W + ^(1)^(0) + ^(0)^(1) + ^(1)^(0) 

= (M2)(i)^(i), (43c) 

= (M2)(i)7y(i). (43d) 

It is easy to recognize that equations ()43a[) and (I43b[) are 
not independent from each other, since at the Alfvenic 

point Ia, F[flF^f = F[f / Fl^f = {UA/Ue){lA). Hence 
at I A, the unknowns ry'-"-' and rj'^^^ can all be ex- 
pressed in terms of which can be arbitrarily cho- 
sen. Once these four quantities are known, the wave 
quantities ^ and rj can be obtained by the expan- 
sion (|42)) at the two grid points immediately astride 
the Alfvenic point. Equation ([38|) is then integrated 
by using a fourth order Runge-Kutta method both in- 
wards to the co ronal base Rp, and outwards to 1 AU (see 
iCranmer fc van Balleg ooiien.200 5). We then rcscale the 
obtained solution such that \ue\ is fixed at Rq for all w 
to be considered. 

5. NUMERICAL SOLUTIONS 

In this section, the solutions to equation ([55]) with a 
number of different angular frequencies uj are presented. 



cj=10"Vad a;=10"*rad s"' cj=10"°rad s 




1 10 100 1 10 100 1 10 100 

HKe) r(R3) r{Ji„) 



Fig. 2. — Radial profiles of the Fourier amplitudes. Results are 
shown for the real (solid curves) and imaginary parts (dashed) for 
three angular frequencies u) = 10~^ (left column), 10~* (middle) 
and 10~^ rad s~^ (right), (a), (d) and (g) Magnetic fluctuation 
given in terms o{ b/Bi, Bi being the background magnetic field 
strength, (b), (e) and (h) Proton velocity fluctuation iip. (c), (f) 
and (i) Alpha velocity fluctuation iia ■ Note that panel (g) uses a 
scale different from (a) and (d). 

All these solutions have the same amplitude for the elec- 
tron velocity fluctuation \Uf.\ = 10^/2 km s~^, or equiv- 

1/2 

alently the time average (ul) ' = 10 km s-^ Note 
that such a value is only about 1/3 of the upper limit 
of the nonthermal vel ocity amplitude deriv ed from line 
width measurements ()Baneriee et al.l [19981) . With this 
choice at the coronal base, we can avoid awkwardly large 
wave amplitudes in interplanetary space, a natural con- 
sequence of the assumption that no wave dissipation is 
present. Nevertheless, the base amplitude can be seen 
as a free scaling parameter for equation ((38|) . The rel- 
ative deviation of the non-WKB from WKB results is 
independent from this choice, even though the absolute 
magnitudes are affected. 

Figure [2] presents the radial profiles of the real 
(solid lines) and imaginary (dashed lines) parts of the 
Fourier amplitudes with three different angular fre- 
quencies, u = 10-3 (left column), 10"'' (middle) and 
10^^ rad s"^ (right). Since the modulus of the magnetic 
fiuctuation |6| spans several orders of magnitude, b/Bi is 
plotted instead of b in panels (a), (d) and (g). Note that 
panel (g) uses a scale different from (a) and (d). Panels 
(b), (e) and (h) depict the proton velocity fluctuation Up, 
while panels (c), (f) and (i) give the alpha one Ua- 

The most prominent feature of Figsl2^, [2)3 and[2j; for 
UJ ~ 10~^ rad is that the radial dependence of the 
real and imaginary parts of fluctuations exhibits a clear 
oscillatory behavior. This is particularly true beyond 
^ 10 Rq. Further inspection of the region r >1Q Rq in- 
dicates that the magnetic and fluid velocity fluctuations 
are well correlated, as would be expected from the eigen- 
relation ^ obtained in the WKB limit. The WKB na- 
ture is further revealed by examining the phase relation, 
i.e., the displacement between the real and imaginary 
parts, for any of the three Fourier amplitudes. For in- 
stance, for the Up profile the nodes in the real part corre- 
spond well to the troughs or crests in the imaginary part. 
Now if examining the envelopes of the Fourier amplitude 
profiles, one can see that \b\/Bi (Figl2^) tends to increase 
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with r while \up\ tends to decrease beyond say 10 Rq. On 
the other hand, |mq| exhibits a non-monotonic behavior 
and possesses a local minimum at ^ 68 Rq. Moreover, 
the magnitude of Ua is considerably smaller than that 
of Up beyond 0.3 AU. This is understandable in light of 
equation ([25|) since the local phase speed of Alfven waves 
Uph is close to the alpha flow speed Ua (see FiglTjD). 

The Fourier amplitudes for lu = 10"'* rad (the mid- 
dle column) also possess an oscillatory behavior, how- 
ever the exact correlation between the magnetic (Figl2jl) 
and ion velocity (Figsl2t and[2f) fluctuations is gone, as 
is the exact phase relation between the real and imag- 
inary parts. Furthermore, one can see that the mag- 
nitude of the alpha velocity fluctuations (Figl^f) in in- 
terplanetary space is substantially larger than that for 
u) = 10~^ rad s~^. Now let us move on to the right 
column for which u) — 10"^ rad s~^. Distinct from the 
preceding two columns, the oscillatory feature disappears 
altogether. Instead, the real and imaginary parts of the 
Fourier amplitudes evolve slowly with radial distance r. 
We note that such a transition from wave-like to quasi- 
steady depe ndence on r with decreasing u; was expl ored 
in detail b y 'Heinemann & Olbert' fl98ff) (see also iLoul 
ll993HMacG rcgor & Charbonncau 1994). That there ex- 
ists a critical frequency below which the transition oc- 
curs was interpreted in terms of the coupling between 
inwardly and outwardly propagating waves. Although in 
the present study, the propagation in opposite directions 
has not been explicitly separated, one can see that for 
a realistic 3-fluid solar wind, there also exists a similar 
critical frequency which may be approximated by 

~ UraAK^TA) (44) 

(cf. Eq.(55) of lHeinemann fc Olbertill980f) . For the cho- 
sen background flow, we find Wc ~ 3.8 x 10~^ rad s^*. 
The numerical experiments have confirmed the validity 
of this approximation. 
Figure [3] shows the radial evolution of the time- 

1 /2 

averages of the magnetic (6^) as well as the fluid veloc- 

1/2 

ity fluctuations (u\) [k — p,a) for the three angular 
frequencies ui — 10^^ (solid lines), 10"** (dashed lines) 
and 10~^ rad s~^ (dash-dotted fines), and for the WKB 
case (dotted lines) as well. The WKB results are eval- 
uated from equations (^5)) and ([50]) by using the same 

1/2 

value for (ttg) at Rq as in the numerical solutions. It 
is obvious that for all the frequencies considered, the pro- 
files demonstrate substantial deviations from the WKB 
one, not only in the absolute values but also in spatial de- 
pendence. A WKB-like spatial dependence is recovered 
only in profiles with oj = 1Q~^ rad s~^ and r > 10 Rq or 
those with u = 10^"* rad s^* in regions r > 120 Rq. 

1/2 

It is interesting to see that, the profiles are 

strikingly similar in the region r < 10 Rq for all three 
frequencies. This behavior can be understood as follows. 
First of all, it can be readily shown from equation (|38p 
that 

(M2 - l)(|ep)' = 2Re(^^n)|^P + 2Re{F,2Vn- 

For the solutions with all three frequencies, it turns out 
that the first term on the RHS always dominates the 
second one in the inner corona. In addition, in the first 
several solar radii, the flow is very sub-Alfvenic, i.e.. 



10-' I ' — ' — — 
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■pi 
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1 10 100 

Fig. 3. — Time-average of the wave-related fluctuations, (a) 
Magnetic fluctuation (h'^'^^^'^\ (b) and (c) Proton and alpha veloc- 
ity fluctuations (u'^'^^^^ and (u^)"'^''^. Numerical results for three 
angular frequencies as well as for the WKB estimates are given in 
different line styles as indicated in panel (a). 

M| < l.Xme < 1, and ^ « -1. It follows that Re(Fii) 
is dominated by -ZilnRBi)' « (lni?B,)'. Therefore |CP 
is approximately proportional to l/{R^Bf). From the 

1 /2 

definition of ^ fEQ. (|39p ). one can see that {by - 1/R. 

That is, {b'^y^^ shows little frequency dependence. On 
the other hand, from equation ([50]) it follows that in the 

WKB limit {b^y^^ - for M| < 1 since Uph « Ua. 
Hence the difference between the profiles for the three 
frequencies and the WKB one in the region r < 1.7 Rq 
reflects the fact that the mass density has a scale height 
less than i?/4. Now consider the portion r > 50 Rq 
where M|. ^ 1. One may see that the slope of the 

1/2 

(6^) profile for lu — 10~^ rad s~* becomes similar 
to that for oj = 10~^ rad s~* (or equivalently the WKB 
one) asymptotically. On the other hand, the profile with 
uj = lO"'^ rad s~^ has a flatter slope than the WKB one. 

1 /2 

This is understandable by noting that the (6^) profile 
for u! = 10~^ rad s~* can be rather accurately repre- 
sented by the zero- frequency solution (|35ap . By noting 
that pUm/Bi is a constant and M|. ~ r^^ for r > 50 Rq, 

~ 1 /2 

one can see that |6| and therefore <^5^) ~ r^*. On the 
other hand, for the WKB profile one can see from equa- 
tion ([301) that asymptotically (6^)^^^ p^/^ r^^/^ 
since Uph ^ Um + UA. 
Moving on to FigslSa and one may notice that 

1/2 

the magnitude of the ion velocity fluctuations ^u^) 
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(fc — p, a) have nearly the same values at i?© as assumed 

1/2 

for (wg) , which can be expected from equation ([7]) to- 

1/2 

gether with the fact that <C i?; for r w i?©. Let 

us examine FigOj in some detail. One can see that be- 
tween 1 i?0 and 1 AU, the WKB estimate always exceeds 



1/2 

the computed (up) value. Furthermore, for r < 5 R, 

1/2 

the (wp) profile for uj = IQ^'^ rad s^^ differs little from 
that for Lu = 10"'' rad s" 
larger than that for w 

5 Rq for instance. One finds that (up)'''^ is 83.9 km s~ 
for the WKB estimate, but 44.6 and 41.9 km s"^ for 
CO = 10~^ and lO"'' rad s~^, respectively. As to the case 

1/2 



■0: 



however they are significantly 
= 10~^ rad s~^. Take values at 

,1/2 



= 10-5 rad s-i, (ul) ' is 17.2 km s 



In other 

words, the WKB estimate yields a value that is 1.88, 
2 and 4.88 times the computed values for lu = 10~^, 
10~'^, and lO"'' rad s^^ respectively. At 1 AU, the dif- 
ference between the WKB estimate and computed results 
is also substantial. One finds for to = 10~^, 10~^, and 
10-5 rad s-\ the WKB estimate is 1.83, 1.18 and 1.32 
times the values obtained numerically. An interesting 

1/2 

feature of the ^u^) profile for r > 0.3 AU is that it 
approaches a constant asymptotically. This can be ex- 
plained by noting that the magnitudes of ion velocity 
fluctuations for uj — 10^^ pad s^^ can also be r oughly 
represented by the zero- frequency solution (j35bp . (The 
zero frequency solutions psp can be accurately repro- 
duced if UJ is further red uced. ) At distances R » Ra 
where Af^ ^ 1, equation (|35bp can be written as 



AnR 



Up + CUa 



UmUp 



1 

Ml 



i?2 



,(45) 



where the constant ^ = [paUa) / {ppUp) is the ion mass 
flux ratio, and is 0.19 for the supposed background solu- 
tion. In the region considered, the first term is found to 
dominate the second in the square parentheses. As a re- 
sult. Up ~ RUap ~ rllA- Since Ua ~ r~^, one finds that 

1/2 

\up\ and therefore {up') approach a constant asymp- 
totically. 

Let us proceed to Figl3t, from which one can see that 

1 /2 

the (u^) profiles also deviate considerably from the 

5i?©,the(u2)'/' 



WKB estimate. In particular, below ' 

1/2 

profiles behave in a fashion similar to (wp) , however 
the difference between the case for lu = 10~^ rad 
and the WKB estimate is even larger. For instance, at 
5 Rq , the WKB estimate is about 10 times the value with 
On the other hand, for 




(Ho) 



Fig. 4. — Time-average of (a) the wave energy density Euj, and 
(b) wave energy flux density .F^ as evaluated from equation l|20|l . 
Numerical results for three angular frequencies as well as for the 
WKB estimates are given in different line styles as indicated in 
panel (a). 



equation (|35b[) that 
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,(46) 



in which the first term in the square parentheses, as in 
the case for (^UpY^^ , is found to dominate the second. 



1/2 

Therefore (mq) ~ Ci'^p)^'^ asymptotically. As a re- 

1 /2 

suit, (u^) should show little spatial dependence. 

Figure [3] presents the radial distribution of (a): the 
wave energy density and (b): the energy flux den- 
sity Fw for the three angular frequencies and WKB ex- 
pectations. In both FigsHK and the proflles for 
LU = lO"'^ rad s"^ approach a WKB-like behavior when 
r > 10 i?Q. When lu = 10'^ rad s'S a WKB- 
like spatial dependence of and F^^ can be seen for 
r > 80 -Rq. Now let us consider Figg^ first. At 1 i?©, 
for LU — 10^5 rad s^^, E^^ is larger than the WKB result, 
but for LU = 1Q~^ or 10~^ rad s~^, it is smaller. This 
behavior is determined by the magnetic component in 
Euj (Eg. ipO)) ) since the particle part is nearly frequency 
independent. Note that in the WKB limit, the wave en- 
ergy is equally distributed between the kinetic and mag- 
netic energies, J2k Pk ("fe) /2 = (^^) /(Stt). Hence at the 



.1/2 



LU = 10~5 rad s~'^. On the other hand, for cu = 10" 
(10^''^) rad s~^, one finds the ratio of the WKB value 
to that numerically derived is 1.89 (2.1). At larger dis- 
tances, one may notice the local minimum in the pro- 
flle with UJ = 10^^ rad s"^ at 0.3 AU, consistent with 
Figl2t. Furthermore, one can see that the values for 
LU — 10^5 rad s^^ become an order of magnitude larger 
than the WKB one. At 1 AU, for lu = 10"^, lO"'', and 
10~5 rad s~^, the values obtained numerically are 0.55, 
1.44 and 11.9 times the WKB one, respectively. In addi- 
tion, similar to (u^)^''^, (m^)^''^ also approaches a con- 
stant value asymptotically. This is understandable since 
in the region R ^ Ra where Af| ^ 1, one finds from yields 



coronal base, for w = 10 and 10 rad s 



the ki- 



netic energy is larger than the magnetic one, whereas for 
UJ — 10~5 rad s"-', the tendency is reversed. It is in- 
teresting to note that among the proflles for the three 
frequencies the highest uj corresponds to a proflle that 
deviates the most from the WKB proflle in the region 
r < 2 Rq. In the region r > 30 Rq, one can see that 
for UJ = 10"'* rad s~^ the E^, profile becomes WKB-like 
asymptotically. However, for uj = 10"^ rad s~^ the slope 
of the proflle is flatter than the WKB one. This 
can be understood since for R ^ Ra, the WKB limit 
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On the other hand. 
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Fig. 5. — Time-average of the wave acceleration exerted on 
(a) protons a„_p and (b) alpha particles a„ ^ as evaluated from 
equation lllll l. Numerical results for three angular frequencies as 
well as for the WKB estimates are given in different line styles as 
indicated in panel (a). In panel (b), the wave acceleration changes 
sign when the profile crosses the horizontal line. 



for R 3> Ra, the zero- frequency result yields that is 
largely determined by the magnetic energy, which leads 
to Ey, ~ ^ r-2. 

Inspection of Fig|3|3 reveals that the the spatial depen- 
dence of the profiles of the energy flux density is re- 
markably similar in the region r < 2 Rq. In contrast, in 
the same region, FigH^ shows that the E^ profiles with 
different u> have rather different r dependence. This be- 
havior of Fu, stems from the fact that the loss of the wave 
energy flux in the form of the work done on the plasma 
is negligible. In other words, is diluted only by the 
flux tube expansion below 2 Rq (cf. Eq. fTS)) '). For in- 
stance, for uj — 10^'^ rad s~^ at 2 Rq, J2k PkUk0.w,k is 
43 erg cm^^ s^^/Rq, amounting to only 3.5% of F^/r. 
As a result, {Fyj/Bi)' « 0, or Fyj cx Bi. That is, in 
this region the spatial profile of Fyj should show little 
frequency dependence. This is found to be true for all 
the frequencies considered, and for the WKB result as 
well. On the other hand, in the region R ^ Ra, one 
can see that for ut = 10~* rad s~^ the F^ profile be- 
comes WKB-like asymptotically, whereas the profile for 
UJ = 10~^ rad s~^ decreases more slowly than the WKB 
one. It turns out that asymptotically for the low fre- 
quency waves Fw is dominated by the Poynting vector 
which evolves like ~ pUmR^M^ ~ M|. ^ r~^. In the 
WKB limit, however, the ratio of the contributions of 
particles to the Poynting vector is roughly 1/2. As a 
result, F„-C/p,,|6|2 ~|&|2^r-3._ 

From the perspective of solar wind modeling, one may 
be more curious about the feedback from the waves to 
the background flow. To this end. Figure \5\ presents 
the radial distribution of the acceleration exerted on (a) 
the protons a^^^p and (b) alpha particles a^^a- One can 



see that for cu ~ 10^^ and 10^** rad s^^, the wave ac- 
celeration is less effective than that in the WKB limit 
throughout the computational domain. However, in the 
low frequency case uj — 10^^ rad s^^, a^^p exceeds the 
WKB expectation considerably beyond 22.8 Rq. Nev- 
ertheless, in all cases the wave force tends to accelerate 
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Fig. 6. — Possible modification of ion flow speeds due to the 
wave acceleration. Although the speeds Up^corr and Uc,,com eval- 
uated from equation I I47I I. are not computed self-consistently, their 
deviation from the background flow parameters represents the sig- 
nificance of the wave acceleration in the force balance for the two 
ion species. Numerical results for three angular frequencies as well 
as for the WKB estimates are given in different line styles as indi- 
cated. 



the protons, i.e., a^.p > everywhere between the coro- 
nal base and 1 AU. However this is not the case when 
the alpha particles are concerned. It can be seen from 
FigO) that the wave acceleration aw,a is negative in the 
interval between 22.7 and 66.1 Rq in the WKB limit as 
well in the case uj = 10"'^ rad s~^. This coincidence 
of the positions where aw,a changes sign stems from the 
fact that beyond ~ 10 Rq, the wave is WKB-like. On 
the other hand, the a^ ^ profile for uj = 10^** rad s^^ 
changes sign at 21.8 and 61.8 Rq. These locations are 
slightly different from their counterparts in the WKB 
limit. When it comes to uj = 10~^ rad s~^, the a^^a 
profile does not show any resemblance to the WKB ex- 
pectation. In particular, the waves start to decelerate 
the alpha particles even in the inner corona: a^^a is 
negative everywhere beyond 3.4 Rq. An interesting as- 
pect of the wave-induced acceleration is that asymptoti- 
cally both aw,p and a^^a approach constant values. This 
is understandable since in the zero-frequency limit, in 
the region where R ^ Ra, both Up and Uk show lit- 
tle radial dependence. It then follows from the expres- 
sions and that for ion species k {k = p,a), 
0"w,k ~ u^(lni?)' — UkUk{\^R)' b/ Bi oc R'uk. Since in 
the region considered, the line of force is nearly perfectly 
radial, one can see that R' is a constant. As a result, 
o,w,k should approach a constant asymptotically. 

The extent to which the wave forces may alter the 
ion flows can be obtained only through a self-consistent 
mo deling by using, e.g., the ite r ative approach adopted 
bv lMacGregor fc Charbonneaul ()1994D : The wave equa- 
tion ((38)) and the solar wind equations dH) to pO)) in- 
corporating the wave contribution are solved alternately 
until a convergence is met. As a first step, however, we 
may simply evaluate the ion speeds Uk,corr {k = p, a) 
corrected for the wave force, i.e.. 



k.corr 



u. 



(47) 



Figure [5] presents the radial distribution of both Up^corr 
and Ua,corr for all the frequencies considered. The back- 
ground flow speed proflles are also plotted for compar- 
ison. One may see that, with the present choice of the 
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Li and Li 



wave amplitude, the waves have neghgible effects on the 
ion acceleration below the Alfven point. Beyond the 
Alfven point, the effects introduced by the waves on 
the speed profiles become more important, especially 
in the low frequency case. As a matter of fact, for 
ui — 10^^ rad s'-^ the corrected proton speed Up^corr 
reaches 785 km s^^ at 1 AU where the background value 
is 648 km s~^. As for the alpha speed, Ua.corr becomes 
negative beyond 100 Rq due to the significant decelera- 
tion exerted on the alphas by the waves. Of course, such 
a situation will not appear in reality. What will hap- 
pen is that the protons are accelerated whereas alphas 
are decelerated by the low frequency waves until the ions 
move at nearly identical speeds. The net effect of low- 
frequency waves is thus to limit the speed difference be- 
tween the protons and alpha particles at large distances. 
However, it should be pointed out that in these regions 
the net work done by the low-frequency wave on the so- 
lar wind as a whole is negligible: J^k PkUkaw,k ~ and 
Fw is nearly divergence free as discussed in reference to 
FigH 

6. SUMMARY AND CONCLUDING REMARKS 

This study has been motivated by the apparent lack 
of a non-WKB analysis of Alfven waves in a multi-fluid 
solar wind with differentially flowing ions. To be more 
specific, this study is concerned with the propagation 
of dissipationless, hydromagnetic (angular frequency lo 
well below ion gyro- frequencies), purely toroidal Alfven 
waves that propagate in a background 3-fluid solar wind 
comprised of electrons, protons and alpha particles. Az- 
imuthal symmetry is assumed throughout. No assump- 
tion has been made that the wavelength is small com- 
pared with the spatial scales at which the background 
flow parameters vary. The wave behavior at a given to 
is governed by equation (j38p . which is derived from the 
general transport equations in the five-moment approxi- 
mation. The Alfvenic point, where the combined Alfven 
Mach number Mt = 1 (cf. Eq. ([Ml) ). is a singular point 
of equation p8|) and a regularity condition has to be 
imposed. For the other boundary condition, we impose 
a velocity amplitude of 10 km s^^ at the coronal base 
(1 Rq). For the given background model of a realistic 
low-latitude fast solar wind, equation ([38| is integrated 
numerically for three representative angular frequencies 
LO = IQ-^, IQ-"^ and 10~^ rad s"^ to yield the radial 
distribution of the wave energy and energy flux densities 
as well as the wave-induced acceleration exerted on ion 
species. 

The first conclusion concerns the applicability of the 
WKB approximation. Between 1 Rq and 1 AU, the 
numerical solutions show substantial deviation from the 
WKB expectations. Even for the relatively high fre- 
quency uj = lO^'^ rad s~^, a WKB-like behavior can 
be seen only in regions where ^ 10 Rq. In the low- 
frequency case llJ = 10~^ rad s~^, the computed profiles 
of wave-related parameters show a spatial dependence 
that is distinct from the WKB one, the deviation be- 
ing particularly pronounced in interplanetary space. In 
the inner corona r < 4 Rq, the computed ion velocity 
fluctuations are considerably smaller than the WKB ex- 
pectations in all cases, as is the computed wave-induced 
acceleration exerted on protons or alpha particles. As for 
the wave energy and energy flux densities, they can be 



enhanced or depleted compared with the WKB results, 
depending on lu. 

The second conclusion is concerned with how the wave 
acceleration may alter the background flow parameters. 
In reference to FiglH it is found that with the current 
choice of base wave amplitude, the wave acceleration has 
little effect on the force balance for protons or alpha par- 
ticles in the corona. That is, one has to invoke processes 
other than the non-WKB wave acceleration to accelerate 
the ions out of the gravitational potential well of the Sun. 
However, at large distances beyond the Alfvenic point, 
low-frequency waves may play an important role in the 
ion dynamics, with the net effect being to equalize the 
speeds of the two ion species considered. 

Strictly speaking, the separation of the flow into fluc- 
tuations and a time-independent background implies 
that the waves are linear. However, one may have 

1/2 

noticed that the wave amplitude (6^) at 1 AU for 
Lu = 10~^ s~^ rad s~^ is substantially larger than 
the background poloidal magnetic field strength (cf. 
FigE^)- That the transverse magnetic field dominates 
the poloidal one demands a careful examination of the 
nonlinear effects other than the wave-induced accelera- 
tion. In particular, one needs to look at the generation 
of secondary waves and structures by the primary Alfven 
waves through the sour ce te r ms in the momentum equa- 
tion. As discussed bv iLoul (jl993l ). in the case of ideal 
MHD these source terms decrease sufficiently fast with 
radial distance asymptotically. Consequently, the first- 
order wave amplitudes are valid provided that the am- 
plitude imposed at the coronal base is sufficiently small. 
The basic picture is expected to be the same even if a sec- 
ond ion species is included, although a similar discussion 
in the 3-fiuid framework will be complicated by the rich- 
ness of wave m odes due to the differe ntial proton-alpha 
streaming fe.g.. iMcKenzie et al.lll993[) . 

As has been mentioned in the introduction, Alfven 
waves are dissipated in some way, and this dissipation of 
the primary waves should be described self-consistently. 
A possibility to do this is to perform a full Elsasser 
analysis extended to the multi-ffuid case, and to express 
the dissipation in terms of the amplitudes of counter- 
propagating waves. We note that this already compli- 
cated issue will become even trickier considering the ne- 
cessity to apportion the dissipated wave energy among 
different species. 

In closing, we note that the low-frequency waves may 
also be important for outfiows from stars other than the 
Sun. For instance, in the radiatively driven stellar winds, 
these waves will provide a further channel of momen- 
tum exchange between passive ions and line-absorbing 
ions in addition to the Coulomb fric t ion. T his possibility 
was first pointed out bv lPizzo "eTalllligsl in connection 
with the effects of stellar rotation. Due to the clear re- 
semblance between the low-frequency Alfven waves and 
stellar rotation (cf. section [221), their discussion also ap- 
plies to the case where the star persistently emits Alfven 
waves with frequencies lower than the critical one de- 
fined by equation Consequently, the mass loss rate 
may be significantly altered. A quantitative study of this 
effect is beyond the scope of the present paper though. 
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